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We study the exchange coupling mediated by itinerant carriers with spin-orbit interaction by both 
analytic and numeric approaches. The mediated exchange coupling is non-collinear and its spatial 
trends depend on the Fermi surface topology of the itinerant carriers. Taking Rashba interaction as 
an example, the exchange coupling is similar to the conventional Ruderman-Kittel-Kasuya-Yosida 
type in weak coupling. On the other hand, in the strong coupling, the spiral interaction domi- 
nates. In addition, inclusion of finite spin relaxation always makes the non-collinear spiral exchange 
interaction dominant. Potential applications of our findings are explained and discussed. 
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I. INTRODUCTION 

The central theme of spintronics is to manipulate 
the ext ra spin de grees of freedom in condensed matter 
systems ^ * 2 * 3 * 4 * 5 * 6 ^ . One of the proposals is to manipulate 
the spin polarization of electrons by the spin-orbit (SO) 
interactional such as the Rashba interaction from the 
structural inversion asymmetry at a surface or interface 
in semiconductors^. To make these proposals feasible in 
spintronics, it is important to achieve a flexible tuning of 
the interaction strength. So far, many surprising works 
have demonstrated how to enhance the strength of the 
Rashba interaction either in a semiconductor or in some 
metal spin-splitting surface states. For instance, by dop- 
ing Bi into GaAs, it was shown that the spin-splittin g gap 
is significantly enhanced in the alloy GaAsi-xBiJ 10 * 11 !, 
which opens up an alternative approach to enhance the 
Rashba coupling in semiconductors. Furthermore, it was 
shown recently that the surface states of Bi/Ag show 
a giant Rashba-type spin splitting due to the surface- 
potential barrier^. In addition, it was also demonstrated 
that, by doping Pb into Bi/Ag, the Fermi energy of the 
surface states is tunable 13 . These exciting discoveries 
open up the possibilities to manipulate the Rashba cou- 
pling over a wide range of interaction strengths. There- 
fore, it is important and exciting to explore the related 
physics properties in both weak and strong regimes where 
the Fermi surface (FS) topologies are different. 

Note that the Rashba interaction in a two-dimensional 
electron gas gives rise to two different FS topologies 
that may play a significant role in dictating the re- 
sponse functions. The classical Dyakonov-Perel (DP) 
mechanism of spin relaxation^ due to Rashba cou- 
pling shows different signatures when the FS topology 
changes. In weak coupling, the DP relaxation mech- 



anism leads to exponential decajEl But, in strong 
coupling (or equivalcntly, the low-density regime), the 
DP relaxation mechanism changes into a peculiar power 
law^Sl. Furthermore, in strong coupling, the cancella- 
tion of the spin-Hal l con ductivity for spin-conserving 
momentum scattering } 16 * 17 ! is not related to the vanish- 
ing of the vertex function, as happened in the weak 
regime. Instead, the cancellation is between the on- 
Fermi surface and off- Fermi surface contributions^. An- 
other example is the disorder-induced localization due to 
the Rashba coupling 19 . For a two-dimensional electron 
gas with strong Rashba coupling, it behaves effectively 
as an one-dimensional system^, which results in infi- 
nite bound states even for a short-range potentia l 19 * 21 ! 
The topological change of the FS also plays an inter- 
esting role in superconductivity^. It was shown that 
the critical temperature of the superconductor can be 
tuned by the Rashba coupling, suggesting materials with 
stronger Rashba coupling are good candidates for en- 
hanced superconductivity 2 ^-. Finally, in our previous 
workP^, we demonstrated that the exchange coupling 
mediated by itinerant carriers with Rashba coupling is 
strongly tied up with the FS topology. We showed nu- 
merically that the non-collinear exchange coupling re- 
veals two spatial trends depending on the FS topologies - 
one is the RKKY-like oscillatory behavior and the other 
is the spiral. 

Here we extend our previous workP^ by a complemen- 
tary analytic approach within the path-integral formal- 
ism. This method is particularly powerful in describing 
the asymptotic behavior of the mediated exchange cou- 
pling and provides a comprehensive understanding of the 
general trends. In the weak Rashba regime, the exchange 
coupling consists of the RKKY oscillations at shorter 
length scale and a spiral background 25 . On the other 



2 




FIG. 1: Noncollinear exchange coupling mediated by itinerant 
carriers with Rashba interaction and finite spin relaxation. 
The hard magnet on the left is pined while the direction of 
the soft magnet on the right is determined by the effective 
exchange coupling. 



hand, the RKKY interaction is suppressed when enter- 
ing the strong coupling regime. It is rather nice that the 
analytic results agree with previous numerics. We also 
use a phenomenological approach to investigate the ef- 
fect of the finite spin relaxation in the itinerant carriers. 
Since the finite relaxation rate breaks the time-reversal 
symmetry, it leads to the suppression of the RKKY inter- 
action. Thus, the non-collinear spiral interaction stands 
out and becomes dominant. 

The rest of the paper is organized as follows. In Sec. II, 
we revisit the ferromagnetic/normal/ferromagnetic tri- 
layer junction mediated by a two-dimensional electron 
gas with Rashba interaction. The band structure and 
the Fermi surface topologies at different carrier concen- 
trations are explained in detail. In Sec. Ill, we derive 
the mediated exchange coupling within linear response 
theory. After integrating out the angular part, the resul- 
tant formula can be implemented in numerics to calcu- 
late the effective exchange coupling. In Sec. IV, we use a 
complimentary analytic approach to derive the effective 
exchange coupling. In Sec. V, the effects of finite spin 
relaxation are discussed in both weak and strong Rashba 
regimes, followed by a brief conclusion. 



II. LINEAR RESPONSE THEORY 

Consider a ferromagnetic/normal/ferromagnetic 
(F/N/F) trilayer magnetic junction (TMJ) as shown 
in Fig. [I] where the intermediate layer contains a two 
dimensional electron gas (2DEG) with the Rashba cou- 
pling. The 2DEG in the intermediate layer is described 
by the HamiltoniarP, 



H R = / d 2 r ^ 
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where is the strength of the Rashba interaction and 
m* is the effective mass of the itinerant carriers. The two- 
component spinors ^\ 'F are the creation/annihilation 
operators for the iti nerant car riers. After integrating out 
the itinerant carrier j 2 3 ! 24 ! 29 ! 30 ^ the exchange coupling be- 



tween the ferromagnetic layers is described by an effective 
Heisenberg Hamiltonian, H e g = £V JijS l L S 3 R . Within 
the linear response theorjPIEl the mediated exchange 
coupling Jjj is prop ortional to the static spin susceptibil- 
ity tensorFllElffl 



Xij(r, 9) . 
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It is rather remarkable that the symmetry arguments 
make the angular dependence explicit and reduce the nu- 
merical task down to evaluation of four real scalar func- 
tions, go(r), gi{r), 32(f) and ho{r). After integrating out 
the angular parts, the first three functions can be casted 
into the following integral forms, 
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where the dispersion for the particle with momentum k 
and chirality A = ±1 is etx = (k — A&r) 2 /(2m*) — Eft. 
Here we introduce the Rashba momentum k R = m*^ji 
and the Rashba energy En = kj 1 /2m*. The oscillatory 
kernels (after integrating out the angular parts) in the 
integral are 
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where J n (x) is denoted the Bessel function of the first 
kind. For the Rashba Hamiltonian, its simplcncss gives 
rise to an extra relation ho{r) — go (r) + 32(f) beyond 
the symmetry argument. As a result, we only need to 
evaluate three independent scalar functions. 

Supposing the ferromagnet on the left of the trilayer 
magnetic junction is aligned along the z-axis, we are in- 
terested in the mediated exchange coupling proportional 
to Xiz( r i 9 — 0), where r is the width of intermediate 
layer. The induced moment is captured by the spiral an- 
gle (shown in Fig. [T]), which can be expressed in terms 
of Xxz and \zz, 
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Therefore, by evaluating g a (r) numerically, we can study 
the spatial trends of the spiral angle 4>x{f). 

In our numerical computations, we will change the pa- 
rameter kfi/kp, where fc# = m*^ R and kp = ^/2m*ep 
to calculate the mediated exchange coupling in different 
regimes. For ballistic carriers, rj — > is assumed in all cal- 
culations. For realistic materials or in the weak Rashba 
regime, we choose the spin splitting A# = 2kp^n= 5 
meV and the Fermi energy ep — 60 meV. Or equivalently, 





FIG. 2: The non-collinear angle <j) x in the weak Reshba regime 
withfcfl/fcF = 0.02 << 1, 77 — * and Lf is defined as jf^. 
The green dotted line is represented the upward tendency 
resulted from the Rashba effect. Notice that the lumps at 
discontinuous points are come from numerical error. 



it corresponds to the Rashba coupling 7^ = 8.91 x 1CP 12 
eVm and the carrier density n^D — 1-25 x 10 12 cnr 2 . 

The results are shown in Fig. [2] In the weak Rashba 
regime k^/kp -C 1, the spatial trends of the angle <p x (r) 
oscillates near the vicinity of ±7r/2 and is very similar 
to the ordinary RKKY interaction. The main difference 
is that, in the presence of the spin-orbit coupling, there 
is an upward background. The oscillatory behavior is 
characterized by Lp = 2n/2kp as in ordinary RKKY 
oscillations. It is important to emphasize that the spa- 
tial trend is drastically different from the spin-precession 
argument in Datta-Das SFETSI because the quantum in- 
terferences from different patches of the Fermi surface 
cannot be ignored. 

The RKKY oscillation with an upward trend can be 
understood in a simple picture. Taking the asymptotic 
limit kpr ^> 1, the reduced spin susceptibility along the 
radial direction Xab{r), where a,b = x,z, can be well 
approximated as ID Ra shba system. Applying a lo- 
cal gauge transformation^ 1 ^, U(r) = e ~ lkRraV / 2 , the 
Rashba Hamiltonian can be mapped into the ID free 
electron gas with the well-known RKKY spin suscepti- 
bility. Since the local gauge transformation is nothing 
but the local rotation about the y-axis with the spiral 
angle <fr(r) = k^r, the reduced susceptibility is approx- 
imately the usual RKKY oscillation twisted by a local 
spiral transformation, 



FIG. 3: The non-collinear angle 
regime with kn/kF = 3.5 and 77 
spiral angle is Lf — n /(kit + Uf). 



in the dilute density 
0. The period of the 
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where the summation over the repeated index c = x, z 
is implied. The gauge argument explains why our nu- 
merical results resemble the RKKY oscillation but with 
a gradual spiral background. 

By increasing the strength of Rashba coupling or low- 
ering the electron density, we enter the strong coupling 



regime k^/kp 3> 1. In this regime, the angle <fix{r) 
ceases to oscillate and rotates gradually with the junc- 
tion width as shown in Fig. [3j At first glance, the spiral 
evolution is similar to the spin-precession argument in 
Datta-Das SFET. But, it is not exactly so. The semi- 
classical spin precession is described by the length scale 
L s = 2n/(2kfj). But, as Fig. [3] shows, the spiraling 
length in our numerics is Lf — 2n/(2kji + 2k p) < L s . 
Only in the extreme strong coupling where kp/ku <C 1, 
the semiclassical picture will be completely correct. 

It is important to mention that the momentum 2fcfj + 
2k p connects patches of the outer Fermi surface which 
are related by time-reversal symmetry. The period of the 
spiral angle indicates that the exchange coupling mainly 
contribute by the interference of the time-reversal sec- 
tion of outside chirality. Since the mediated exchange 
coupling is spiral- like, many phenomena, in this regime, 
related to the RKKY exchange coupling should be recon- 
sidered. For example, a tunable non-local spin control 
in a coupled quantum dot system via a gate voltage has 
been beautifully dem onstrate d in the 2D electron gas sys- 
tem of semiconductor^ 27 * 28 ^. The authors showed that 
the Kondo effect is suppressed by a nonlocal RKKY-like 
interaction which can be used to control the quantum 
dot spin. If the mediated exchange coupling is no longer 
RKKY-like, it is interesting to see how the non-collinear 
coupling will reshape the idea. 



III. PATH INTEGRAL FORMALISM 

In this section, we derive the effective exchange cou- 
pling analytically. Taking the asymptotic limit, the spa- 
tial trends of the non-collinear exchange coupling become 
rather clear. Within the path integral formalism, the ef- 
fective coupling is obtained by integrating out the itin- 
erant carriers and retaining the action to the quadratic 
order. After some algebra (details can be found in Ap- 
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pendix A), the effective Hamiltonian can be written as a 
contour integra l 25 * 34 * 35 * 



H. 



eff 
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where /i is the chemical potential, R = Ri — R2 and 
r\ represents an infinitesimal imaginary energy, and tr 
means a trace over the spin degrees of freedom. In 2D 
Rashba system, the retarded Green's function in momen- 
tum space takes the form, 
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Fourier transforming it back to the coordinate space^^, 
the retarded Green's function along the x-direction is 



G(R; e + i0+) = G (R; e)l + G X {R- t)a y 
with Go and G\ defined as 
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where R = |R| and q = \Jk 2 + k 2 R with k\ — 2m* e. 
Hn^ is the Hankel function of the first kind, defined as 
Hn\z) = J n (z) + iY n (z) with J n and Y n being the Bessel 
functions of the first kind and second kind. Notice that 
the Bessel or Hankel functions have a branch cut along 
the negative real axis and the infinitesimal i0 + is impor- 
tant in the spiral regime. We substitute the retarded 
Green's function into Eq. (10 1 and obtain the effective 



Hamiltonian along the x-axis, 

H cS = ASi ■ S 2 + B(Sx x S 2 ) y + CSfSl 
where coefficients A, B, and G are 
-J 2 
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The hard magnet is aligned along the z-axis, S = Sfz, 
so that we can rewrite the effective Hamiltonian as 



U cS = ASfSi + BSfS%. The function A(R) represents 
the collinear exchange coupling and B(R) for the non- 
collincar part. By computing these two functions, the 
spiral angle can be determined. 

Let us start with the weak Rashba coupling. In this 
regime, we have q > k R , and by ignoring the slight contri- 
bution of k e < 0, the Green's functions can be simplified 
to be 
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With the asymptotic form of the Hankel functions, 
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we end up with 
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The effective Hamiltonian is thus greatly simplified in 
the asymptotic limit, 

Tifi^^F^R) [cos(2k R R)S[SI + sia(2k R R)SfS^],(2A) 
where the range function^ 



F 1 (R) 



J 2 m* sin(2q F R) 
IStt 2 R 2 ' 



(25) 



Note that the range function is identically the same as 
the usual two-dimensional RKKY range functio n^, ex- 
cept that kp is replaced by qp = y/2m*ep + k 2 R . There- 
fore, as clearly demonstrated in Eq. (24), the mediated 
exchange coupling is similar to the usual RKKY with a 
non-collinear spiral background as obtained numerically 
in Fig. [2J 

Now we turn to the strong coupling regime. In this 
regime, because (q — k R )R < 0, we need to be cau- 
tious about the analytic properties of the Hankel function 
Hn\z). On the complex plane, Hn\z) has a branch cut 
on the negative real axis. Thus, we need to pay extra at- 
tention to the i0 + factor. Note that the Hankel function 
of the first kind along the negative real axis can be con- 
nected to the Hankel function of the second kind, 



H^(-x + i0 + ) = (-l) n+1 H^(x + tO + ), 



(26) 



1 2 ) 

with the definition H\ (z) = J n (z) — iY n (z). By this 



connection, we can rewrite the Eq. (13 1 and alAn as 
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FIG. 4: The non-collinear angle (f) x in the critical regime with 
hii/ftF = 1 and r\ — > 0. The spiral length at the critical point 
is 7r/2fc_R = L 3 /2, only half of that in the strong coupling 
regime. 



In the asymptotic limit, the second kind of Hankel func- 
tions are 



the numerical results are shown in Fig. [4] It is clear 
that the trend belongs to the spiral type. However, pay 
special attention to the length scale of the spiral. In 
the strong coupling regime, the spiral length is roughly 
Lj. Right at the critical point (k R — kp), the spiral 
period is n/(2k R ) — L s /2, which follows the rule in the 
strong coupling regime. This factor of two can be under- 
stood as the more general formula for the length scale 
L+ = ir/(k R + hp) in the strong coupling limit and at 
the critical point. The analytic calculations can be done 
by substituting q = k F = k R into Eq. (13 1 and Eq. (14 1. 



The asymptotic forms of the Green functions near the 
Fermi surface arc 



Go — G\ 
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2^irk R R 
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(35) 



and the effective Hamiltonian is 



n c cS oc [cos(4k R R)S^SI + sin{4k R R)SfS$] , (36) 

with the spiral period 2ir/4k R = L s /2, as we mentioned 
above. 
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Simple algebra leads to the following results, 
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Finally, the effective Hamiltonian takes the form, 
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with the range function 
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where Si(x) — J* m f^dt is the sine integral function and 
is positive for x > 0. In the asymptotic limit, Si(x) 
reaches the constant value it/2 with minor oscillations. 
Thus, the spatial trend of the mediated exchange cou- 
pling is dominated by the spiraling part with minor os- 
cillations. This is drastically different from the results 
in the weak Rashba regime. Again, our analytic results 
agree with those shown in Fig. [3] obtained by numerical 
methods. 

The same calculation can be done for the critical case 
between the weak and the strong coupling regimes where 
one of the Fermi surface disappears. For comparison, 



IV. FINITE SPIN RELAXATION 

So far, we have revealed the interesting transitions be- 
tween weak, strong and critical regimes assuming the me- 
diating carriers are perfectly ballistic. In realistic mate- 
rials, the spin relaxation rate is finite due to many mech- 
anisms. Thus, it is interesting to study how the spatial 
trends of the mediated exchange coupling changes when 
the carriers are no longer perfectly ballistic. To fully 
address this issue, one needs to rederive the effective ex- 
change coupling using a modified single-particle propa- 
gator. While this is certainly an interesting direction 
for future studies, we peek into the problem by intro- 
ducing a phenomenological spin relaxation ration within 
the linear response theory. It is motivated by the obser- 
vation that the inclusion of finite spin relaxation breaks 
time-reversal symmetry and shall suppress the RKKY in- 
teraction. Therefore, we expect the non-collinear spiral 
interaction should be enhanced and becomes dominant. 

For convenience, we introduce the parameter k, q = 
y/2m*r]. In Fig. [5] we choose the parameters k v /kF = 1 
and k R jkp = 0.02, and we can see that the RKKY 
pattern is indeed suppressed significantly due to the fi- 
nite relaxation rate. Furthermore, as we move closer 
to the critical regime ( but still in the RKKY regime) 
with k R /kp — 0.2, even a moderate spin relaxation 
k v /kF = 0.5 will almost wash out the signature of RKKY 
oscillations and the spiral trend becomes rather robust, 
as shown in Fig [6] Therefore, it seems that the finite 
spin relaxation will suppress the RKKY oscillation and 
enhance the spiral interaction. 

In conclusions, we apply both analytic and numeric 
methods to study the exchange coupling mediated by the 
itinerant carriers with Rashba interaction. We demon- 
strate that the Fermi surface topology greatly alters the 



6 



APPENDIX A: 




20 R(nm) 40 



FIG. 5: Suppression of RKKY oscillations in the presence of 
spin relaxation, the RKKY phase with kn/kF = 0.02 and 
k v /k F = 1. 




R(nm) 



FIG. 6: Suppression of RKKY oscillations in the presence 
of spin relaxation with closer to the phase transition with 
kn/kF = 0.2 and k v /kF = 0.5. 



property of mediated exchange coupling in 2D Rashba 
gas system. When the Rashba energy is much smaller 
than the Fermi energy, the mediated exchange coupling 
is almost the RKKY oscillation plus a slight upward 
trend. On the other hand, if the Rashba energy is larger 
than the Fermi energy, the mediated exchange coupling 
shows a spiral pattern. Apparently, the Fermi surface 
topology plays a crucial part about the mediated ex- 
change coupling. We suggest that, in the quantum-dot 
systems! 251221221 with dilute carrier densities, the spiral- 
like exchange coupling should be considered instead of 
the usual RKKY interaction. 
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from the National Center for Theoretical Sciences in Tai- 



DERIVATIONS FOR EFFECTIVE 
HAMILTONIAN 



Within the path integral approach, we investigate the 
effective coupling between two ferromagnets S\ and S 2 
mediated by itinerant carriers, as illustrated in Figjl] 
Consider the Zcncr model, but only included two im- 
purity spin, the Hamiltonian can be represented as 



H — H + Hi, 



(Al) 



where 



H =f d D r J d D r'J2^(r)h (ry)M>(r'), (A2) 
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JSiir) ■ - 
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Here we ignore the correction effects of the itinerant car- 
riers, and represent the unperturbed Hamiltonian Hq as 
two fermionic operators. ^ is a two-component spinor for 
the itinerant carriers and Sj(r) = SjS(r — Rj) denotes 
the impurity spin at the position Rj in the D-dimension, 
and J represents the coupling between the itinerant and 
the localized spin densities. 

In order to write down the partition function in path 
integral form, we introduce the Grassmann variables 
4>, (j> for the itinerant carriers^ and the coherent state 
for the impurity spins^, which is denoted as S\Cl) = 
SCl\Cl), where spin orientation is denoted here by £1 = 
(sin 9 cos ip, sin 9 sin ip, cos 9) and the length of impurity 
spin vector by S. In path integral language, the parti- 
tion function is represented as 



X>[fti,fi 2 ] / V 



■ f" dTL(fii,n 2 
Jo 
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where the integral denotes the summation of infinite time 
section. The Lagrangian can be writte as 

+H(n 1 ,n 2 ;4>,cb), (A5) 
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where the first and second terms are the dynamics of the 
impurity spins and itinerant fermions respectively. Inte- 
grate out the itinerant carrier and arrive at an effective 
description for localized spin, 



Z = / V[n u n 2 ] e --s e ff(n 1 ,a2) j 
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with the effective action 



Scs= [ dr [ d D r V (0 7 |<9 T |n z ) - In [dct^r 1 ] .(A7) 
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Here we use the notation, 



the effective action after Fourier transformation is 



g- 1 ^!,^) = Go 1 + 5g-\n 1 ,n 2 ), (as) 



Qq 1 = d T l - h . 



(A9) 
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By utilizing the relation, ln[det M] = tr[ln M], the second 
term of the effective action can be expanded, 

In [detg- 1 ] = tr [in (G^ 1 + SG' 1 )] , 

= tr [ln^ 1 ] + tr [in (l + Go^G' 1 )] , 

= tr [inGv 1 ] - t r J2^ (-GoSG- 1 )" ( An ) 

n=l 

The zeroth order term does not contain any spin opera- 
tors and can be neglected when computing the effective 
spin-spin interaction. The first-oder term is linear in the 
spin operator and vanishes since the itinerant carriers 
consider here are not polarized, 

tr^S- 1 ] =~ JdTjd D v Go(r,r;v,T + 0+) 

x tr |<T • fi/(r,r)| = 0, (A12) 



7=1,2 

where r + + is for the creation operator. 

The second-order term contains two spin operators and 
is relevant for the effective exchange coupling, 

S' cS = itr [GoSG-'GoSG- 1 ] 



4 



dr / dr' 
[i Jo 



X 



'trj^-O^r)] Go(Rut;R 2 ,t') 
?-f1 2 (r')]0o(R2,T , ;Ri,T)}. (A13) 



By truncating the expansion to second order, the ef- 
fective action for the two impurity spins is S c g = 

Jodr J d D rJ2i =h2 (^i\ d rPi) + S' eS . The first term of 
this effective action, known as Berry phase, describes the 
dynamics of the impurity spins. Here we assume that the 
dynamics of the itinerant spins is much faster so that we 
can ignore the retardation effect between two impurity 
spins. In other words, we take the stationary limit of the 
impurity spins £l/(r) = Si and ignore the Berry-phase 
term. 

Now we would like to write down the effective action 
in the Fourier space of the imaginary time^, 



-ik n (t— T ) 



Go{ik n ) (A14) 



where the summation is over all Matsubara frequencies 
k n = — — for fermions. Making use of the identity 



/3 



P 



dr e l(fe '"- fc " )r = 6„ 



(A15) 



J 



3 ■ S± 
3 ■ S 2 



g (Ki,K2;ik n ) 
5o(R2,Ri; **«)}. (A16) 



Notice that by neglecting the retardation effect, the poles 
of the Green's function, Go{z), on ly show up on the real 
axis. This property let Eq. ( A16 1 to be rewritten in a 



retarded Green's function form. 

First, we can demonstrate the stationary RKKY 
effect between two impurity spins mediated by free 
fermions from Eq. (AI6l. By Fourier transforming 



Eq. (AI6) to the momentum space and fixing the di- 



rection of one impurity spin to the z-axis, we have 

2SfS 2 , and thus the equation 



tr 



(3 ■ Six) [3 ■ S 2 ) 



can be rewritten as 



S e fi — 







dr H eS , 



(A17) 



with H eS = J(R i2 ) SfSi and 



J 



j (ft \ _ " / d q f d k p .qRn 



^g (k + q;zfc„)g (k;ifc n ). (A18) 







ik,, 



Here we use the translation invariant property of free car- 
riers and assume no spin polarization of this system. The 
summation of Matsubara frequencies for fermions can be 
computed and represented as the particle-hole propaga- 
tion. After some calculation, we end up with 



J(Ri 2 ) 



(2nY 



jq R12 



x(q), (A19) 



where the particle-hole propagation is denoted as 



x(q) = / £ 



(2tt) d 



n H£k+q) - M£k) 



- ^k+q + "1 



(A20) 



where ^ is the eigenenergy of free fermions and np(^) 
is the Fermi-Dirac distribution function. As shown as 
the Eq. ( |A19 1 , we obtain the mediated RKKY exchange 
coupling strength between two impurity spins. 

Next, we will replace the imaginary Green's functions 
in Eq. (A16) by the retarded ones. The connection be- 



tween these two types of Green's functions can be s een 
in the complex planed The kernel of the Eq. (A16l in 
the complex plane, 



/(*) =tr{ 



3 ■ Si 



g (Ri,R 2 ;z) 



3 ■ S 2 



0o(R2,Ri;*)}, (A21) 
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FIG. 7: Contours in the complex plane. The black points 



denotes the Matsubara frequencies z 



■ (2n + l)7T 





for fermions. 



has poles on the real axis. If we consider a contour inte- 
gral, 



I = 



dz 



2m 



,f{z) n F (z) e z0+ , 



(A22) 



where the contour path is shown in the Fig|7] it includes 
two parts, the upper and lower of the complex plane. The 
contour integral can be rewritten, 



j = Res 

z — ik n 



f(z)n F (z)e z0+ ] = --3^2 f (ikn). (A23) 



Here the Fermi-Dirac distribution is represent as n F (z) — 
(eP z + and has poles at z = i(2n + l)n/(3, n — 

1, 0, 1, • • • , where the residues at these poles are 

Res [n F {z)] = lim Z -=^- = -\. (A24) 



integral can be simplified to two line integrals very close 
to the real axis, 



,00+ir, dz 

I = lim i I - — : n F (z)f(z) 



lim 



— oo+irj 



2m 



-°°~ tri dz 
2wi 



n F (z)f(z)\, (A25) 



OO — 17] 



+0+ ./_oo 271"? 



n F {u) [/(w + *T ? )-/(w-tT?)].(A26) 



Combining Eq. ( A23 1 and Eq. ( A26 1 , we obtain the useful 
identity, 



V/(ifc n ) = - / dr / 



2tti 



n F (uS) 



[f(w + if!) - f(LJ - ir,)} , (A27) 

where the first term of the right hand side is represented 
by the retarded Green's function and the second term 
is denoted the advanced Green's function. By using the 
identity, 

tr [ABCD]* = tr [B t A t £) t C t ] = tr [A t Z> t C t B t ] , 

(A28) 

and the relation between advanced and retarded Green's 
function, [Go(R-2, R-i; w + iff)]* = Go(Ri, R2; u> — iv)i we 
can obtain the second term as complex conjugate of the 
first. After reshaping the effective action, we end up with 
the effective Hamiltonian of the Eq. (Al7l, 



H e ff = 



Because in the |R| — > 00 limit, the factor n F (z)e z0+ ap- 
proaches zero whether Re(z) > or not, the contour modified effective Hamiltonian comes to be the Eq. (10 1 



-J f°° ( \( ^ \ 

— — J doj n F (oj) Im | tr ya-Sij 

G (Ri, R 2 ; w + iri) (ff ■ S 2 ) G (R 2 , Ri; u + irj)\ } (A29) 

At T = 0, the Fermi function is a step function. We can 
simply replace n F [ui) by 1 and change the upper limit 
of the integral to the chemical potential fi. Furthermore, 
if the unperturbed system is translational invariant, the 
effective Hamiltonian is just the function of the relative 
distance and we denote R42 = Ri — R 2 - In the end, the 
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